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Abstract 

Recently, a strongly enhanced xy anisotropy of magnetic excitations was observed in YBa2Cu30 y 
(YBCO y ) with y = 6.45 and T c = 35 K [Science 319, 597 (2008)]. Unlike the observation in 
YBCOg.6 an d YBCOg.85) the anisotropy grows to be pronounced at lower temperature and at 
lower energy, and is not suppressed by the onset of superconductivity. We propose that the effect 
of singlet pairing is substantially reduced in YBCO6.45. This reduction concomitantly enhances 
an order competing with singlet pairing, a strong tendency of the so-called d-wave Pomeranchuk 
instability, leading to the magnetic excitations observed experimentally. 

PACS numbers: 74.25.Ha, 74.72.Bk, 74.20.Mn, 71.10.Fd 
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The spin-charge stripe order was widely discussed in high-T c cuprates. [1|] In this scenario, 
the two-dimensional Cu02 plane is assumed to have the instability of one-dimensional charge 
order, the so-called stripe pattern, in the antiferromagnetic (AF) background with a 7r 
phase shift across the charge stripe. However, the charge order signal was observed only for 
limited materials with specific hole doping rates and had a broad spectrum with very weak 
intensity, implying that the charge order is usually not well developed. Kivelson, 

Fradkin, and Emery introduced a new concept of electronic nematic order, [6j] which was 
envisaged as the melting of charge stripe order, [7] where the translational symmetry is 
recovered but the orientational symmetry is still broken. 

On the other hand, an alternative root to yield electronic nematic order was found in 
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minimal models of high-T c cuprates such as the t-J s|, [lo| and Hubbard [11 
models. These models show a tendency toward rf-wave type Fermi surface deformations 
(dFSD), the so-called c?-wave Pomeranchuk instability. The Fermi surface (FS) expands 



along the k x direction and shrinks along the k y direction, or vice versa. The <iFSD state 
has the same symmetry as the electronic nematic state. But the physical origin is different 
from the stripe physics, [6] since the c£FSD is generated by forward scattering processes of 
quasiparticles, not by fluctuating charge stripes. 

The rfFSD competes with superconductivity. [if]] In the slave-boson mean-field^ and 
variational Monte Carlo flOj] analyses of the t-J model, superconductivity becomes dominant. 
Yet the system still has sizable ciFSD correlations. [171] producing a giant response to a 
small external anisotropy, e.g., due to a lattice structure and anisotropic strain, si, 10] This 
theoretical insight yields a promising scenario 18j to understand the pronounced anisotropy 
of magnetic excitations observed in untwinned YBa 2 Cu 3 0, / (YBCO y ) with y = 6.5, 19] 



6.6,[20, 



211 ] and 6.85, 20j where the lattice yields a small xy anisotropy. 



Quite recently, a very strong anisotropy was observed in magnetic excitations in un- 
twinned YBCO-6. 45.(22|] The observed spectra were qualitatively different from typical ob- 
servations in Y-based cuprates. (i) The anisotropy starts to increase below ~ 150 K and 
saturates below ~ 50 K, in contrast to the case of YBCO-6.6 20, |2lJ where the anisotropy 
is most pronounced at relatively high temperature and is reduced at low temperature, (ii) 
The magnetic excitation spectra are hardly affected by the onset of superconductivity and 
in fact a gap feature is not found in the magnetic excitations. This is sharply different from 
the well-known observation in Y-based cuprates, where a broad spectrum around q = (n, n) 
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becomes very sharp and is strongly enhanced below the onset temperature of superconduc- 
tivity, leading to a resonance at relatively high energy, followed by suppression of low energy 
spectral weight. 
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24j These peculiar magnetic excitations in YBC0 6 . 45 are not captured 
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implying failure of the 



even qualitatively by simple application of the theory in Ref. 
standard RVB mean-field theory in a strongly underdoped region. 

Why do the qualitative features of magnetic excitations change so drastic in YBCO6.45? 
Since YBCO6.45 lies closer to the AF instability, strong AF fluctuations might play an 
important role. However, even for the lower hole-doped system YBCO6.353, and thus much 
closer to the AF instability, the experimental data 25( showed that magnetic correlation 
length is still short ranged. In this paper, we show an integral role of <fFSD correlations, 
which are concomitantly enhanced by the suppression of singlet pairing, leading to results 



very similar to the experimental observation. 



22] 



We analyze the bilayer t-J model on a square lattice, a minimal model for Y-based 
cuprates, [26] 

H = - 2J trC^Cr'a + ^ J r S r • S r > (1) 

r,r',o- ( r ! r '} 

defined in the Fock space with no doubly occupied sites. The operator cJ CT (c rcr ) creates 
(annihilates) an electron with spin a on site r, and S r is the spin operator. J T (> 0) is a su- 
perexchange coupling between the nearest neighbor sites along each direction, r — x,y, z. We 
take into account hopping amplitudes t T between r and r' up to third-nearest neighbors, and 
the direction of r' — r is represented by r. The Hamiltonian (CQ) is analyzed in the standard 
slave-boson method in the t-J model by introducing the so-called resonating-valence-bond 
(RVB) mean fields. However, in contrast to the standard RVB theory, we assume the mean 
field of singlet pairing A to be zero, in order to mimic possibly substantial suppression of 
singlet pairing in YBCO6.45 as we will discuss later in detail. The formalism and the band 



parameters are otherwise the same as those in Ref. 
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where we successfully discussed the 
pronounced anisotropy of magnetic excitations in YBC06.85 [20] and YBC06.6- 20, 21] We 
introduce 5% xy anisotropy into the hopping integrals t T and 10% into the superexchange 
coupling J T , twice as large, as imposed by the superexchange mechanism. The anisotropy 
then can be parameterized by a single parameter a = —0.05. The minus sign indicates that 
the band parameters are more enhanced along the y direction, parallel to the direction of 



the CuO chains, which indeed reproduces 
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28] qualitatively the same anisotropic dis- 
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FIG. 1: T dependence of the anisotropy of renormalized band, 2Ai/Aio, for several choices of 5 
for a = —0.05. 



tribution of magnetic excitations as observed previously. 20j, |2lj Interestingly the negative 



sign of a is a 



excitations. 29] 



so implied in a spin spiral state to understand the anisotropy of magnetic 



The most important quantity is the renormalized in-plane nearest-neighbor hopping t T , 
which is given by i T = t T 5 + |J r Xr with r = x,y; 5 is the hole density and Xr is the RVB 
bond order along the r direction. The order parameter of the <iFSD is, then, defined as 
Xd = {Xx — Xy)/2- Figure [1] shows the temperature (T) dependence of the anisotropy of 
renormalized band 2Ai/i , where At = \i x — t y \/2 = \(t x — t y )5/2 + \{J X — Jy)Xd\ and 

= (tx + t y )/2. The anisotropy is strongly enhanced at lower T, especially for lower 



5. This enhancement comes 



rpm the development of Xd, namely the underlying <iFSD 
flp| The presence of the anisotropy in t T and J T , however, 



instability in the t-J model, 
smears the singularity associated with the (iFSD instability, which would appear below 
T w 0.03J(0.05J) for 5 = 0.08(0.07) and disappear for 5 > 0.09 if the input anisotropy is 
zero. 

The irreducible dynamical magnetic susceptibility Xo(q, oj) reads 

. . 1 v — > tanh ^fe — tanh !^i, q 

= + » + P) 

where N is the total number of (bilayer) lattice sites and T is a positive infinitesimal. We 



calculate magnetic excitations in the renormalized random phase approximation (RPA) 30, 

1 + ^(q)Xo(q, w) 

Here 

J (q) = 2r(J x cos q x + J y cos q y ) + J z cos q z (4) 
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FIG. 2: (Color online) Left-hand panels: q maps of Inrx(q, uj) for a sequence of T in 0.6tt < g^, q y < 
IAtt for <5 = 0.08, uj = 0.05J, and a = —0.05. Right-hand panels: corresponding results for the 
bare anisotropy, namely without dFSD correlations. 



with a renormalization factor r. While r = 1 in the plain RPA, the value of r is reduced in 
the renormalized RPA such that the AF instability appears in a doping region comparable 
with actual materials. We set r = 0.4, which confines the AF instability to 5 < 0.068. 

We calculate x(q, uj) numerically by choosing V = 0.05 J. This choice of V is mainly due to 
numerical convenience, but simulates damping of electrons by static defects in real materials 
and broadening due to limited energy resolution in neutron scattering experiments. Since 
neutron scattering measurements for untwinned YBCO are performed for the odd channel, 
we focus on this channel (q z = tt). Considering that the hole density in YBCO6.45 is estimated 
as 8.5%, [22] we fix S = 0.08. 

The left-hand panels of Fig. [2] show two-dimensional q maps of the imaginary part of 
^(q, uj) at uj = 0.05 J for a sequence of temperatures. At T = 0.01J, significant spectral 
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FIG. 3: (Color online) (q, uj) map of Imx(q, u) at T = 0.01J for J = 0.08 and a = -0.05; the q 
directions are shown in the inset. 

weight is centered around q = (tt, 7r) and is strongly elongated along the q x direction. 
This pronounced anisotropy originates from the strong enhancement of the original band 
anisotropy due to <iFSD correlations (Fig. [1]). The anisotropy of magnetic excitations is 
reduced with increasing T. In particular, the distribution of the spectral weight becomes 
nearly symmetric around (n, 7t) at T = 0.30 J. To clarify the <iFSD effect, we show the 
corresponding result for the bare anisotropy effect by switching off dFSD correlations in the 
right-hand panels of Fig. [2j where we impose the same anisotropy a = —0.05, but assume 
Xd = 0. Although the value of At does not depend on T, we see that the anisotropy of 
Inr\;(q, a;) depends weakly on T and becomes visible with decreasing T. But the obtained 
anisotropy [Figs. E^a') andH^b')] is much weaker than the results due to <iFSD correlations 
[Figs. [2](a) and [2(b)] . Hence the underlying rfFSD correlations are particularly important to 
produce a strongly enhanced anisotropy at lower T. 

To see the u dependence of the anisotropy, the (q, u) map of Imx(q, u) is shown in Fig. [3] 
at low T. While the anisotropic distribution is still discernible for a high energy where the 
spectral intensity is substantially reduced, a pronounced anisotropy appears especially in a 
low energy region. 

The evolution of the spectral intensity of Imx(q, uj) is shown in Fig. H]as a function of uj 
for several choices of T at q = Q = (tt, tt). At high T, Imx(Q, uj) shows a broad maximum at 
a moderate energy. With decreasing T, low-energy spectral weight is substantially enhanced, 
yielding a peak structure at low energy. 

When g?FSD correlations strongly enhance a small xy anisotropy of the original band 
dispersion, a strong anisotropy of IniY(q, uj) is induced with a characteristic T dependence. 
In YBCO6.45, the experimental data 22j suggest that fermions responsible for magnetic exci- 
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FIG. 4: uj dependence of Imx(Q, oS) for several choices of T for 5 = 0.08 and a = —0.05. 

tations feel a rather small mean field of singlet pairing in the sense that neither a gap feature 
nor a clear impact of the onset of superconductivity is observed in neutron scattering. This 
special experimental situation is mimicked by setting A = in the slave-boson scheme of 
the t-J model. We have found that thisphenomenological treatment well captures the most 



salient features observed in YBCO 



6.45- 



221 ] (i) the strongly enhanced anisotropy of Im%(q, uj) 



at lower T for low lo (Figs. [2] and [3]) and (ii) the enhanced spectral weight of Imx(Q,a;) at 
low uj for low T (Fig. H]). 

n 

On the other hand, in the standard slave-boson mena-field theory of the t-J model[18], 
the onset temperature of singlet pairing is higher than that of the rfFSD, and the singlet 
formation suppresses the susceptibility of the <iFSD, which then does not diverge. [8| The sus- 



ceptibility, however, remains appreciable 17] so that a small external anisotropy introduced 
in the t-J model is strongly enhanced even in the singlet pairing state. [s] This effect is more 
pronounced for a lower doping rate. Such an enhanced anisotropy nonetheless turns out 
not to be sufficient to capture the anisotropy of Imx(q, to) observed in YBC06.45[22] even at 



qualitative level for a realistic choice o 
at relatively high T (Fig. 20 in Ref. 



parameters. In particular, the anisotropy is enhanced 



18l ). Moreover the experimental observation 22J that 



the onset of superconductivity hardly affects magnetic excitations is difficult to be captured 
in existing microscopic theories for cuprates. 

Inclusion of singlet pairing is not so straightforward, and something beyond the existing 
theories happens in YBCO6.45. Here focusing on the context of the present paper we dis- 
cuss several possibilities that (iFSD correlations still remain strong enough to be compatible 
with the present result even in a full calculation. First, since the g?FSD instability con- 
cerns discrete symmetry breaking, the spontaneous symmetry breaking is allowed at finite 
temperature even in the exact analysis of the two-dimensional model, in contrast to super- 
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conductivity and antiferromagnetism. Hence suppression of the ciFSD due to fluctuations is 
less crucial than the latter two, in favor of strong correlations of the (iFSD. In the presence 
of a small external xy anisotropy, then, the anisotropy can be strongly enhanced as we have 
seen in Fig. [TJ The feedback of the enhanced anisotropy appears as further reduction of the 
magnitude of singlet pairing compared with the mean-field value in the slave-boson theory 
of the t-J model. Second, YBCO6.45 has an orthorhombic crystal structure, but is close to 
the continuous transition to a tetragonal structure. The ciFSD order parameter couples to a 
phonon related to the orthorhombic-tetragonal structural transition, which contributes to an 
enhancement of (iFSD correlations. In this case, the phonon spectrum may also show strong 
anisotropy in the orthorhombic phase. Last, there might be a phase segregation into hole- 
poor AF domains with nano-scale correlation lengths and hole-rich superconducting islands 
in YBCO6.45, and the neutron scattering signals originate mainly from the AF domains. Such 
finite-sized domains are often discussed in a strongly underdoped superconducting region in 



yuSR measurements. [32 



While the present result well captures the neutron scattering data in YBC06.45,|22| a 
detailed comparison reveals disagreement about several aspects, which may be mainly due 
to our simple calculation ignoring singlet pairing, (i) A weakly incommensurate structure 
was observed along the q x direction, whereas we have obtained a single peak at (tt, 7r) as 
read off from Fig. [2j In a fermiology scenario, the incommensurate structure in Y-based 
cuprates is explained by the suppression of the commensurate peak due to the development 



of o?-wave singlet pairing. 
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34J Hence the observed weak incommensurate peak can be due 



to the development of relatively weak singlet pairing, (ii) The observed data at uj = 50 meV 
was interpreted to be isotropic. [22J In the present calculation, the anisotropy becomes less 
pronounced for high uj, but is still discernible (Fig. [3]). On the other hand, if we include 

n 

singlet pairing in the slave-boson mean- field theory, it is known [18] that the anisotropy of 
Imx(q, uj) at high uj depends strongly on energy. More detailed experiments as well as more 
complete calculations are necessary to discuss the anisotropy at high uj. 

The present phenomenological calculation suggests that d-wave singlet pairing is substan- 
tially reduced in YBCO6.45, which then concomitantly enhances strong dFSD correlations, 
leading to a magnetic excitation spectrum very similar to the experimental data. {2^ Since 
the magnitude of ci-wave singlet pairing is often discussed to become larger with decreasing 
doping rate, the present study implies an interesting direction to understand the evolution 
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of (i-wave singlet pairing as well as its connection with (£FSD corre 



ations in underdoped 



cuprates. Recent angle-resolved photoemission spectroscopy data[35| is in fact suggestive 
of relatively small <i-wave singlet pairing. It demonstrated two distinct gaps, the amplitude 
of one of which increases for lower doping rates while the other, which shows ci-wave sym- 
metry, decreases. In addition, strong c£FSD fluctuations substantially reduce the lifetime of 
quasiparticles in the antinodal region of the FS while not in the nodal direction as shown 



in Ref. 



3d Hence strong rfFSD fluctuations may play an integral role to understand the 



pseudogap in the strongly underdoped regime of YBCO. Strong rfFSD fluctuations are also 
expected in a scenario of a quantum phase transition into the <iFSD state in a underdoped 



region. 



371 ] Such a scenario is, however, based on the assumption of the existence of robust 



ci-wave pairing, which is different from the present theory where we have invoked the sub- 
stantial suppression of (i-wave pairing due to the competition ly, [17J with the (iFSD in a 
strongly underdoped region. 
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